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ABSTRACT 

The radiation-driven warping instability discovered by Pringle holds considerable promise as the mechanism 
responsible for producing warped, precessing accretion disks in X-ray binaries. This instability is an inherently 
global mode of the disk, thereby avoiding the difficulties with earlier models for the precession. Here we follow 
up on earlier work to study the linear behavior of the instability in the specific context of a binary system. We 
treat the influence of the companion as an orbit-averaged quadrupole torque on the disk. The presence of this 
external torque allows the existence of solutions in which the direction of precession of the warp is retrograde 
with respect to disk rotation, in addition to the prograde solutions that exist in the absence of external torques. 

Subject heading: accretion disks — instabilities — stars: individual (Her X-l, SS 433) — X-rays: stars 


1. INTRODUCTION 

For a quarter of a century, evidence has been accumulating 
for the existence of warped, precessing disks in X-ray binary 
systems. The discovery of a 35 day period in the X-ray flux 
from Her X-I (Tananbaum et al. 1972) was interpreted almost 
immediately as the result of periodic obscuration by a pro- 
cessing accretion disk that is tilted with respect to the binary 
plane (Katz 1973). Katz proposed that the precession was 
forced by the torque from the companion star, but left unex- 
plained the origin of the disk's misalignment. An alternative 
possibility is the “slaved disk" model of Roberts (1974), in 
which it is actually the companion star that is misaligned and 
precessing; the accretion disk (fed by the companion) will track 
the motion of the companion, provided that the residence time 
in the disk is sufficiently short. 

Dramatic evidence for a warped, precessing disk in an X- 
ray binary was provided by the discovery of the relativistic 
precessing jets in SS 433 (Margon 1984, and references 
therein). The systematic velocity variations of the optical jet 
emission, the radio jet morphology, and optical photometry of 
the system all indicate a precession period for the disk of 164 
days; this must be a global mode, as both the inner disk (to 
explain the jets) and the outer disk (to explain the photometry) 
must precess at the same rate. The systematic variation of the 
X-ray pulse profile of Her X- 1 has been interpreted as the result 
of precession of the inner edge of an accretion disk, which 
argues for a global mode in this object also (Triimper et al. 
1986; Petterson, Rothschild, & Gruber 1991). A crucial point 
is that in both Her X-l and SS 433, the direction of precession 
of the warp has been inferred to be retrograde with respect to 
the direction of rotation of the disk (e.g., Gerend & Boynton 
1976, Her X-l; Leibowitz 1984 and Brinkmann, Kawai, & 
Matsuoka 1989, SS 433). 

A number of other X-ray binaries, of both high and low 
mass, show r evidence for long period variations that may in- 
dicate the presence of precessing inclined disks (Priedhorsky 
& Holt 1987, and references therein; Cowley et al. 1991 ; Smale 
& Lochner 1992; White, Nagase, & Parmar 1995): LMC X-4 
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(30.5 days), Cyg X-l (294 days), XB 1820-303 (175 days), 
LMC X-3 (198 or 99 days), and Cyg X-2 (77 days). Given the 
small number of X-ray binaries for which adequate data exist 
to test for the existence of such periodicities, it is evident that 
precessing inclined accretion disks may be common in X-ray 
binaries. 

Although theoretical attempts to understand the mechanism 
responsible for producing precessing, tilted, or warped accre- 
tion disks date from the discovery of Her X-l, no generally 
accepted model has emerged. The original models suggested 
for Her X-l suffer from serious fiaw^s, although both predict 
retrograde precession. The model of Katz (1973) imposes a 
tilted disk as a boundary condition, but provides no mechanism 
for producing this tilt. The slaved disk model (Roberts 1974) 
requires that the companion star rotation axis be misaligned 
with respect to the binary plane; however, the axial till is ex- 
pected to decay by tidal damping on a timescale shorter than 
the circularization time (Chevalier 1976). Other suggested 
mechanisms suffer from the difficulty of communicating a sin- 
gle precession frequency through a fluid, differentially rotating 
disk (e.g., Maloney & Begelman 1997). 

Recently, however, a natural mechanism for producing 
warped accretion disks has been discovered. Motivated by work 
by Petterson (1977) and Iping & Petterson (1990), Pringle 
(1996) showed that centrally illuminated accretion disks are 
unstable to w T arping because of the pressure of reradiated ra- 
diation, which, for a nonplanar disk, is nonaxisymmetric and 
therefore exerts a torque. Further w T ork on Pringle's instability 
w r as done by Maloney, Begelman, & Pringle (1996), who ob- 
tained exact solutions to the linearized twist equations, by Ma- 
loney, Begelman, & Nowak (1997, hereafter MBN), who gen- 
eralized the earlier work to consider nonisothermal disks (see 
below), and by Pringle (1997), who examined the nonlinear 
evolution. Radiation-driven warping is an inherently global 
mechanism; the disk twists itself up in such a way that the 
precession rate is the same at each radius. 

Previous work on Pringle's instability assumed no external 
torques. However, in X-ray binary systems, the torque exerted 
on the accretion disk by the companion star must dominate al 
large radii. In the present Letter we examine the behavior of 
the instability when an external torque is included. In § 2 we 
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derive the modified twist equation and solve it numerically, 
and in § 3 we discuss the solutions and the implications for 
X-ray binaries. 


so that the equation becomes 


d 2 W DW 

x—r + V-ix)—-- 
dx~ dx 


5 {d + c5 0 .r 1 ) IV, 


( 3 ) 


2. THE TWIST EQUATION AND SOLUTIONS 

We use Cartesian coordinates, with the Z-axis normal to the 
plane of the binary system; hence, Z = 0 is the orbital plane. 
As in MBN, we assume that the disk viscosity y, = 

(R/R^y, where R {) is an arbitrary fiducial radius. For a steady 
state disk far from the boundaries, this implies that the disk 
surface density U oc R L In the a- viscosity prescription, <5 = 
2 corresponds to an isothermal disk. We distinguish between 
the usual azimuthal shear viscosity and the vertical viscosity 
v 2 that acts on out-of-plane motions (see, e.g., Papaloizou & 
Pringle 1983; Pringle 1992); the ratio v 2 /v x = 17 is assumed to 
be constant, but not necessarily unity. Assuming an accretion- 
fueled radiation source, we transform to the radius variable 
x = (2 ],2 elr}) (R/R s ) u2 , where e = L/Mc 2 is the radiative effi- 
ciency, and R s is the Schwarzschild radius. The linearized equa- 
tion governing the disk tilt (including radiation torque) for a 
normal mode with time dependence e w ' is then 


D 2 W DW . , „ , 
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equation (10) of MBN. The tilt is described by the function 
W = &e h \ where j 8 is the local tilt of the disk axis with respect 
to the Z-axis, and 7 defines the azimuth of the line of nodes. 
We have nondimensionalized the eigenfrequcncy 0 by defining 
a = 2r) 2 R 2 s (j/€ 4 v {i ; in general, a is complex, with a real part 
a , and an imaginary part a r Negative values of d, correspond 
to growing modes. We also define R 0 such that .v(7? 0 ) — 1 ; 
therefore, R X} = \ ( ri/e) 2 R s . 

Equation (1) assumes that the only torques are produced by 
viscous forces and radiation pressure. In this case, the preces- 
sion of the warping modes must be prograde, i.e., in the same 
direction as disk rotation (MBN). In a binary system, however, 
the companion star also exerts a force on the accretion disk. 
Since the orbital period is much shorter than the viscous times- 
cale, we average the resulting torque over azimuth and keep 
only the leading (quadrupole) terms (e.g., Katz et al. 1982). 
This orbit-averaged torque term will cause a ring of the disk 
that is tilted with respect to the binary plane to process, in a 
retrograde sense, about the normal to the binary plane. (In- 
cluding the time dependence of the torque would result in 
nutation as well as precession, a complication that we do not 
consider here.) Relative to the angular momentum density of 
a ring, the torque term is 
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where R is the distance from the compact object of mass M x , 
and r is the separation between the companion star of mass M 
and the compact object. Writing the ^-dependent coefficient of 
this term as (RfR 0 ) V L ^ 1C quadrupole torque then contributes 
a term k5 0 .v 3 W to the right-hand side of the tw ist equation (I), 


where the quadrupole precession frequency c3 {) has been non- 
dimensionalized in the same manner as the eigenfrequcncy in 
equation (1). 

Equation (3) must be solved numerically. However, an entire 
class of solutions can be derived from the results of MBN, 
which do not include an external torque. Comparison of the 
twist equation with and without quadrupole torque shows that 
equation (3) with a — 0 is formally identical to equation (1) 
(in which <£ {) = 0), if we replace 5 with 8' = 6 — In other 
words, the purely precessing (i.e., real a, since co 0 is real) modes 
with no external torque have the same shapes as nonprccessing 
modes with quadrupole torque, with oj 0 = d„ except that the 
latter correspond to a larger value of the surface density index 
6. The fact that these modes arc nonprccessing, i.e., the warp 
shape is fixed in an arbitrary inertial frame, is rather remarkable, 
as it requires that the radiation torque (w hich attempts to make 
the w’arp process in a prograde sense) and the quadrupole torque 
(which attempts to make the disk process in the retrograde 
direction) precisely balance at each radius. 

We impose the outer boundary condition that the disk must 
cross the Z = 0 plane at some radius, and we impose the usual 
no-torque inner boundary condition. This choice of Z = 0 for 
the outer boundary condition is not strictly correct; what w*e 
take as the disk boundary here corresponds to the circulari- 
zation radius in a real X-ray binary, and the actual outer 
boundary w ill typically be at R iWl ~ 3/? dri .. The proper boundary 
condition then constrains the gradient RW ] DW/DR at /? dlv (J. 
Pringle 1997, private communication). We have modeled the 
outer disk (/? circ < R < R oul \ assuming that - 0 for R> 
R iix and W* = 0 at R altr and have examined how r these solutions 
couple to the inner disk solutions. We find that the gradient is 
always steep at /?. jrc , so that the outer and inner disk solutions 
match up at radii that differ by from the radius of the 

zero, and the lilt declines rapidly to zero for R > /? dlv . Thus, 
the true solutions differ little from the zero-crossing eigen- 
functions calculated here, and we can ignore any minor dif- 
ferences for the purposes of this Letter. A full discussion of 
this important point is given in MBN. 

As in MBN, we separate equation (3) into real and imaginary 
parts and solve as an initial-value problem, iterating to find the 
zero-crossing eigenfunctions. In the absence of any external 
torque, there is a marked change in the behavior of the eigen- 
functions across d — I (see MBN), so we consider here the 
tw'o cases 5 = 0.75 (corresponding to the usual gas pres- 
sure-supported Shakura-Sunyaev disk) and 6 — 1.25. 

In Figure I we plot the location x {) at which the eigenfunction 
returns to the Z = 0 plane (in the terminology of MBN, these 
are the first-order zeros, closest to the origin) and the normal- 
ized precession rate 5,, as a function of the dimensionless quad- 
rupole frequency ca () , for 6 = 0.75. In Figure 2 the same quan- 
tities are plotted for 5 = 1.25. The series of curves are for 
different grow th rates, from a value of o, close to the maximum 
growth rate, down to small growth rates for which the curves 
essentially coincide with the steady state <r, = 0 curve, which 
is not shown. The dotted portions of the curves mark the pro- 
grade modes, the solid portions show the retrograde modes. 
Although there are significant differences in behavior for the 
two values of 6, there are the following similarities overall: 
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Pic j . i. — (a) Top: Location of the disk outer boundary x i} as a function of 
the dimensionless quadrupole precession frequency co ( , for surface-density in- 
dex <5 -■ 0.75. Curves arc for increasing values of the growth rale ~o t from 
right to left; the plotted curves are for ~o t = 0.001, 0.01, 0.02, 0.03, 0.04, 
and 0.05. The solutions with retrograde precession are marked with solid lines, 
the progradc solutions with dotted lines. The maximum growth rate plotted 
(0.05) is very close to the maximum allowed growth rale. ( b ) Bottom: Di- 
mensionless precession rale a, for the same models as in (ri). 


1. There is a maximum value of above which the torque 
from the companion is too large to allow warped modes to 
exist. 

2. For a given value of ca 0 , the solutions are generally double- 
valued, with one prograde and one retrograde mode or two 
prograde modes; the retrograde modes return to the plane at 
larger radius than the prograde modes. 

3. The precession rates of the modes a , are much larger than 
the associated values of a> 0 ; this is unsurprising, since the quad- 
rupole torque fixes the boundary condition at the disk's outer 
edge, which is always at R » R {) (recall x(R {) ) = 1). 

For any single value of the growth rate a,, the disk boundary 
x {) is at a fixed radius for a given value of o> 0 (and choice of 
solution in the double-valued regime). The value of w 0 , in turn, 
is set by the masses and separation of the components of the 
binary. In general this value of .v 0 will not match the actual 
outer boundary of the disk. This implies that in real disks, the 
location of the outer boundary will determine the growth rate, 
provided that the outer boundary falls within the range of radii 
occupied by the warp modes, either prograde or retrograde. As 
can be seen from Figures \a and 2a, the fraction of the (o> 0 , 
.Y u )-plane occupied by retrograde modes is nonncgligible, al- 
though it is smaller than for prograde modes. The retrograde 
modes essentially always have their outer boundaries at larger 
radius than the prograde modes. The precession rates cr, are 
usually comparable for the prograde and retrograde solutions, 
although for 6 = 1.25, |cr,| for the retrograde modes may be 
several times larger than for the prograde modes. (The abrupt 
termination of the curves for these first-order eigenfunctions 


Fig. 2. — («) Top: As in Fig. \a, for 5 — 1.25. From right to left, the plotted 
curves are for a, = 0.01, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.0. 
The maximum growth rate plotted (1.0) is very close to the maximum allowed 
growth rale, (b) Bottom: Dimensionless precession rate a, for the same models 
as in (a). 


is genuine; the higher order eigenfunctions extend to larger a\ 
generally with different values of a r and cr, for a given value 
of w 0 .) 


3. DISCUSSION 

Where do real X-ray binary systems lie in this parameter 
space? In order to quantify our results, we must make a 
specific assumption about the viscosity; we also evaluate the 
results for <5 = 0.75, but there is in fact very little depen- 
dence on <5. The precession rates can be expressed in terms 
of the viscous timescale 7\ tsc ~ RiV K ~ 2/? : /3p, at the critical 
radius R Ki [the minimum radius for instability: typically 
R vr ~ \0 t> 7] 2 {0Afe) 2 M/M Lz) cm; see MBN1, with .v c , = .v(/? cr ): 


7-v,r 6 

o, = — 7T7— o,. 


12 T_(/? tr ) " 


(4) 


Using the a-p rescript ion for viscosity, defining the precession 
timescale to be r r ,. cc = lir/o,, and normalizing to typical values, 
we find that 


25 


M/M 


(e/0. 1) 3 a/0 


— ( HlR \ 2 
A \0.0lA, \0. 1 / 


XT da y s - ( 5 > 


Thus, the expected precession timescales for the disks are from 
weeks to months. 





L46 


MALONEY & BEGELMAN 


Vol. 491 


The quadrupole precession frequency is given by 


oj () = 


3 /CM)' 2 /R,Y n (M\ f 

4 \ r 3 / U UJ 


= 4.9 x 



( 6 ) 


where r u = r/10" cm. Since oi () is normalized in the same 
manner as a, we find that 


w 0 ~2 x 10 


r; 5 lMV r -(HIR\ 2 

£^(a/0.l)/-, , 1 Wy Vo.01 


(7) 


From Figure 1 it is evident that the relevant values of ci t) fall 
in the range in which warped disk solutions exist, with either 
prograde or retrograde precession. Furthermore, the outer 
boundary radius ,v 0 that characterizes these solutions is typically 
R iml ~ 10 5 /? s , which is also the expected value for real X-ray 
binary sytems. 

The estimates of r pfCC and o> 0 are sensitive to the value of the 
radiative efficiency e — extremely so, in the latter case. How- 
ever, we do not expect this to be true in real X-ray binary 
systems. The steep dependence on e in the linear theory esti- 
mate comes from the scaling of the fiducial radius R 0 . Physi- 
cally, however, the important quantity will be the value of oj 
at the disk boundary, and not its value at R u . since it is the 
outer boundary condition that determines the importance of the 
torque from the companion. Thus, wc expect that in reality the 
timescales and frequencies will be less sensitive to e than in 
the above estimates. 

The linear theory of disk warping leaves a number of ques- 
tions unanswered. Although wc have demonstrated that ret- 
rograde as well as prograde solutions exist when a quadrupole 
torque is included, there does not appear to be any reason for 
choosing one mode over another. For a fixed value of u> 0 , the 
retrograde solutions (when both exist) occur for larger values 
of the boundary radius, but this difference is not large, being 
only a factor of ~5-IO at most (see Figs. 1 and 2). In Figure 
3 (Plate LI) we show the shapes of the most rapidly rotating 
prograde and retrograde modes for <5 = 1 .25, for both the steady 
state and fastest-growing modes. (The 5 = 0.75 modes are very 
similar; this difference from the solutions with no external 
torque [MBN1, in which the shapes of the growing modes are 


very different for S> 1 and 5 < 1, is a reflection of the im- 
portance of the quadrupole torque, which always dominates at 
large radius.) The fast-growing prograde modes arc more 
“wound up” (i.e., the azimuth of the line of nodes rotates 
through a larger angle between the origin and the disk bound- 
ary) than the fast-growing retrograde modes, which may result 
in differences in the effects of self-shadowing on the disk. An 
important observational task, given our earlier argument that 
warped processing disks are probably common in X-ray bi- 
naries, w ill be to determine if retrograde precession is the norm 
in such systems or only occurs in some fraction of them. As 
we noted earlier, both of the original suggestions for producing 
processing disks in X-ray binaries require that the sense of 
precession be retrograde. 

We also cannot establish from linear theory that steady long- 
lived solutions actually exist. Work on the nonlinear evolution 
of radiation-warped disks (with no external torques) by Pringle 
(1997) suggests that chaotic behavior can result as a conse- 
quence of the feedback between disk shadowing and the growth 
of the warp. Calculation of the nonlinear evolution of disks in 
X-ray binaries wall be necessary to determine whether steady 
solutions do exist in this case, and if so, under what conditions. 
We have also ignored the effects of winds, which could also 
drive warping (Schandl Sc Meyer 1994; Pringle 1996). How- 
ever, in the case of X-ray-heated wands, this seems unlikely 
to be important, since at the base of the wind, the gas pressure 
is much less than the radiation pressure (e.g., Begelman et al. 
1983). 

Pringle’s instability appears to be very promising as a so- 
lution to this long-standing problem. It is expected to be ge- 
ncrically important in accretion disks around compact objects 
and is inherently a global mode of the disk, so that the warp 
shape prccesses with a single pattern speed. In addition, it 
provides a natural mechanism for producing nonplanar disks 
in X-ray binary systems, so that the torque from the companion 
star (which affects only the out-of-plane portion of the disk) 
is able to influence the disk shape, allowing retrograde as well 
as progradc modes to exist. 
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Fig. 3. — Surface plots showing the shapes of several of the w’arped disk modes for 5 - 1.25. In all cases the amplitude of the warp has been fixed at 20%, 
and the solutions have been plotted to the disk boundary. Top left: The most rapidly rotating (largest a,) steady slate (cr, = 0) prograde mode. Top right: The most 
rapidly rotating steady state retrograde mode. Bottom left: The fastest-growing (maximum -cr,), most rapidly rotating prograde mode. Bottom right The fastest- 
growing, most rapidly rotating retrograde mode. 
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